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We investigate the effect of dynamic electron correlation on high-harmonic generation in helium atoms using 
intense visible light (A = 390nm). Two complementary approaches are used which account for correlation 
in an approximate manner: time-dependent density-functional theory and a single-active-electron model. For 
intensities / ~ 10 14 W/cm 2 , the theories are in remarkably good agreement for the dynamic polarization and 
harmonic spectrum. This is attributed to a low-frequency collective mode together with a high-frequency single- 
electron response due to the nuclear singularity, both of which dominate electron correlation effects. A time- 
frequency analysis is used to study the timing and emission spectrum of attosecond bursts of light. For short 
pulses, we find a secondary maximum below the classical cut-off. The imprint of the carrier-envelope phase, for 
the time-integrated spectral density appears at frequencies above the high-frequency drop-off, consistent with 
previous studies in the infrared A ~ 800nm. 

PACS numbers: 31.15.ee, 32.80.Rm, 42.65.Ky 



I. INTRODUCTION 



The development of sources of coherent ultrashort high-frequency light is a topic of great interest at present, with the prospect 
of illuminating ultrafast electron dynamics 01 El HI]- Over the past decade, high-order harmonic generation (HHG) has been used 
as an extreme ultra-violet light source for a wide range of applications, and is a practical method for producing attosecond bursts 
from intense infrared lasers USUI- In this scheme, an intense laser is focused into a gas target and the high-order nonlinearity 
of the response of the atoms provides a narrow, forward-beam of high-order harmonics. The efficiency of this process depends 
upon the individual atom susceptibility and the phase matching within the gas sample. Experiments based on this technique have 
led to generation of harmonics from rare-gas atoms with an infrared pump laser in excess of the 300th harmonic, using very 
short laser pulses (< 20 femtoseconds) and peak intensities in excess of I ~ 10 14 W/cm 2 [0, 01 . 

Interest is focusing on how these attosecond bursts might be controlled and optimized for applications in ultrafast spectroscopy. 
Availability of intense femtosecond lasers has driven research into the highly nonlinear response of atoms and molecules to 
extreme dynamic fields. The experimental tool of choice has been the TkSapphire laser operating in the near infrared. Recently 
attention has turned to the use of frequency-doubled sources. Naturally, at very high intensities, I > 10 15 W/cm 2 , double 
ionization is a feature of the interaction. However, at intermediate intensities, ionization is a much weaker channel and, to a 
good approximation, the process is governed by single-electron excitation, possibly leading to single ionization. However, even 
in this case, one would expect the electron pair to retain aspects of correlation intrinsic to the unperturbed helium atom. The 
effect of correlation is a fundamental aspect of atomic structure, but its role in terms of hyperpolarizability and high-harmonic 
generation is not fully understood. 

Dynamic correlation in helium can be treated exactly through the direct approach of solving the two-electron time-dependent 
Schrodinger equation (TDSE) 

i^-Mri,r 2 ,t) =JJ¥(ri,r 2 ,t), (1) 
at 

where H is the Hamiltonian operator of the system. The solution embodies the dynamic polarizability that is the source of the 
high-frequency harmonics of interest. This partial-differential equation is intractable (in computational terms) for more than 
two electrons and even in the simple two-electron case, it requires several hundred processors to obtain a numerical solution 
when infrared wavelengths are considered 09D. The reason for these large computational demands is due to the treatment of the 
electron-electron correlation term. In the case of stationary or metastable states of the helium atom, this correlation term can 
be accounted for, in an essentially exact manner, using basis function expansions. However, under transient perturbations, the 
dynamic correlation of the atom is more problematic. Nonetheless, the problem is well known and consequently this equation has 
received a great deal of attention, not least because of its importance in ion-atom collisions in gases and plasmas 01 [TIGS El. 

In this paper, we investigate the role of electron correlation, by using two semi-correlated approaches. More precisely, we use 
two different models which account for correlation in an approximate manner. The methods employed to solve this equation 
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are a time-dependent density functional theory (TDDFT) approach, recently developed to study laser-molecule and laser-cluster 
interactions lfl4ll and a single-active-electron (SAE) model [15]. In this paper, we study the dependence of harmonic generation 
from helium on the intensity and carrier-envelope phase (CEP), both of which have been identified as being important physical 
parameters in the determination of HHG emission properties. 

In Sec. II below, we introduce the TDDFT approach within the exchange-only limit, and provide numerical details for solving 
the corresponding time-dependent Kohn-Sham (TDKS) equations. In Sec. Ill, we describe the SAE method employed. In 
Sec. IV, we briefly mention how the power spectrum is calculated and illustrate how a time-frequency analysis is performed. 
Finally we discuss and summarize the results for HHG in helium. Atomic units are used throughout unless otherwise stated. 



II. TDDFT METHOD 



The TDDFT method has been applied extensively to the study of atomic and molecular systems driven by external laser 
pulses HQ [DIBS USES]. Indeed TDDFT provides one of the most detailed, practical and feasible ab initio approaches 
for tackling many-body problems. The time-dependent formulation of ground-state density functional theory (DFT) B221 12311 was 
provided by Runge and Gross l24ll . who showed that the response of a system of interacting electrons could be obtained from 
that of a set of fictitious non-interacting particles exposed to a time-dependent local effective potential. As well as providing 
a method for studying the time-dependent evolution of an electronic system, it also allows for the calculation of excited state 
properties of static systems. Like TDDFT, many-body effects are in principle included exactly through an exchange-correlation 
functional; in practice the form of this functional is unknown and at best it can only be approximated. 

Our implementation of TDDFT, as applied to a general spin-polarized system of N electrons is set out in reference Hill . In 
the following we describe how our approach is applied to the two-electron atom. The Kohn-Sham wavefunction is written as a 
single determinant of one-particle Kohn-Sham orbitals. Since helium, initially in its 1 S e ground state, is spin degenerate then 
only one Kohn-Sham orbital, ipKs{r, t), is required. The electron spin densities are then equal, i.e. 

n T (r, t) = ni(r,t) = \i) KS (r,t)\ 2 , (2) 

where J., f denotes the spin state of each electron, and so the total electron density is 

n{r,t) =n r (r,t)+n l (r,t) = 2\*p KS (r,t)\ 2 . (3) 

The time evolution of this orbital is governed by the TDKS equation 

d 

«^V'Ks(r, t) = H KS ipKs{r, t), (4) 

where 

H KS (r,t) = -l-V 2 + V ext (r,t)+ [ dr'P^ + V xc (r,t). (5) 
2 J \r — r \ 

In equation (0 the external potential is given by 

F ex t(r,*)=Kon(r)+Kaser(r,t), (6) 

where Vi on (r, t) and V[ asel -(r, t) represent the Coulomb and laser-interaction terms respectively. We consider a linearly polarized 
laser pulse, make the dipole approximation and consider both the length and velocity form of the interaction. In a length-gauge 
description, the interaction term is given by 

V<£(r,t)=r-eE(t), (7) 
where e is the polarization direction. In the velocity gauge, the interaction term is 



vZ}(r,t) = iA(t)e-V, (8) 



where A[t) is the vector potential defined by 



A(t) = A f(t)cos{u; L t + ip), 



(9) 
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and where ojl is the frequency, tp the carrier-envelope phase (CEP) and f(t) the pulse envelope given by 

r sin 2 (nt/T) < t < T, 

/(*) = \ 

\ otherwise 

for a pulse of duration T. With this form of the vector potential, the electric field is 

E(t) = E f(t) sm{u L t + if) - — cos(w L i + tp). (10) 

ujl ot 

1 /2 

The electric field amplitude (E = ojlAq) is related to the cycle-average intensity / by E = (8ttI/c) 1 , where c is the speed 
of light. Similarly we define the ponderomotive energy, Up = E^/(Auj\). 

The third term on the right-hand side of equation <j5j represents the Hartree potential and the fourth term incorporates the 
remaining exchange and correlation effects. This exchange-correlation potential is itself spin-degenerate for helium and can be 
written as 



V xc (r,t) = V XC(r (r,t)= dE * c[n ^ ] 



5n„ 



(ID 



where E xc [n^,ni\ is the exchange-correlation action. 



A crucial element in our model is the functional form of the exchange-correlation potential. While many sophisticated approx- 
imations to this potential have been developed 12511 . the simplest is the adiabatic local density approximation in the exchange-only 
limit (xLDA). The exchange energy functional is then given by 



1 /3 

E x [n h n l ] = -U^-) 5> 4 / 3 (r,i), (12) 



from which the exchange-only potential 



- 6 \ 1/3 

V xa (r,t) = -[-) ?iy 3 M), (13) 



is obtained. While this approximate functional is simple to implement it does suffer from the drawback of containing long-range 
self-interaction errors: the asymptotic form of the potential is exponential instead of Coulombic. The anomalous long-range 
form of the self-interaction potential means that the spectrum of single-particle highly-excited states are incorrect. It follows 
that physical multiphoton resonances, or alternatively (at extremely high intensities) the tunneling process, leading to ionization, 
will be poorly represented. Nevertheless, as we will see, the gross electrical properties of the atom are fairly well reproduced. 

Since we consider a linearly-polarized laser pulse within the dipole approximation, rotational symmetry around the z-axis is 
preserved at all times, and so it is appropriate to solve the TDKS equation using cylindrical coordinates. The electron position 
vector is then given by 



r = p cos tf>i + psinipj + zk. (14) 

Precise numerical details of how the code is implemented are given in ITill . As in ll26ll a finite difference treatment of the 
z-coordinate and a Lagrange mesh treatment of the p-coordinate based upon Laguerre polynomials is employed. The time- 
dependent Kohn-Sham equation is discretized in space using these grid techniques and the resulting computer code optimized 
to run on massively-parallel processors. Several parameters in the code affect the accuracy of the method and these are adjusted 
pragmatically until convergence is obtained. Specifically, these parameters are the number of points in the finite difference grid 
(N z ), the finite difference grid spacing (Az), the number of Lagrange-Laguerre mesh points (N p ), the scaling parameter of the 
Lagrange-Laguerre mesh (h p ), the order of the time propagator (N t ) and the time spacing (At) (3. In all the calculations 
presented here, converged results were obtained using the following parameters: N z = 4485, Az = 0.02, N p = 60, h p = 
0.2027685, Nt = 18 and At = 0.01. Incidentally, these parameters, without any further adjustment, give converged (better than 
1%) TDDFT orbital energies for the complete first-period of elements of the periodic table. 
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III. SAE METHOD 

The single-active-electron (SAE) model 

HOUSED] 

is, perhaps, the simplest and most appealing approach for multiphoton 
ionization in which a single valence electron is released. The model provides results that are cheaply produced, often to a very 
satisfactory degree of agreement with experiment and/or highly -expensive complex many-body calculations. As such, the model 
provides a useful benchmark in the absence of more sophisticated calculations or accurate measurements. We choose to employ 
spherical coordinates with the polar axis along the direction of linear polarization. The radius, r, and polar angle, 9, are treated 
explicitly, whereas </> is treated analytically. Thus the field-free Hamiltonian, Ho, for the single-active-electron, of the helium 
atom initially in its ground-state can be written as 



Ho(r,e) = -\ 



1 d ( o d 



r 2 dr V dr 



1 d ( . „ d 
sin i 



V(r), (15) 



r 2 sin6 dO \ 89 

where the model potential V(r) we chose has the form of that calculated by Tong and Lin lll5ll 

„, v Z + ai e-^ r + a 3 re- a * r + a 5 e- a ° r 

V{r) = , (16) 

r 

in which, Z = 1.0, a\ = 1.231, a,2 = 0.662, a 3 = —1.325, 04 = 1.236, 05 = —0.231 and ae = 0.480. In this expression, 
the long-range monopole Z = 1 is supplemented by short-range corrections, {a;} expressing static correlation. The form of the 
potential was obtained from fitting to a self-interaction-free density functional. The Hamiltonian for the interacting system can 
be written in the form 

H SAE (r, 6, t) = H (r, 9) + VU(r, 9, t). (17) 

The velocity gauge formulation is computationally attractive since the number of partial waves required for convergence is 
greatly reduced compared to the length-gauge. This also provides us with a check on the length-gauge results. 

The TDSE defined by equation (Q]) therefore takes the form 

H SAE {r, 9, *)V>SAE(r, 9, t) = i^ SAE (r, 9, t). (18) 

at 

The SAE wavefunction, V>sae(V, 8, t), is expanded in a direct product of radial and angular functions. The solid-angle normalized 
Legendre polynomials Y®(9) are efficient for the angular dimension, leading to a sparse interaction matrix. The radial coordinate 
is discretized using £?-spline functions l32il . so that 

VsaeM,*) = VVc„,(t)^y»(9), (19) 



r 



with k the order of the splines, and where the expansion is truncated by a maximum angular momentum Z max and the limiting 
number of spline functions n max . This allows us to apply Dirac's variation of constants method to the evolution. Taking inner 
products over the spatial basis functions we obtain a large set of ordinary, coupled first-order differential equations for the 
constants, c n i(t), in the form 

Sc(i) = -iH(t)c(i), (20) 
where the overlap matrix S is due to the non-orthogonality of the B-splines 

S„>i>,nl = 6 VI x / ^ drB k nl (r)B*(r), \n' - n\ < k. (21) 

The interaction matrix is of the form 







/>7T /"^max 

H nn , tnl =27t ddsmO drrB k ,(r)Y l ?(9)H SAE r- 1 B^r)Y l °(9). (22) 

Jo JO 



5 




FIG. 1: (Color online) Generation of high-order harmonics in a helium atom exposed to a five cycle (6.4fs) linearly polarized laser pulse of 
wavelength A = 390nm for two different models: a TDDFT model (black lines) and a SAE model (red lines). The plots on the left-hand-side 
correspond to a laser intensity of I = 1 x 10 14 W/cm 2 while the plots on the right-hand-side correspond to a laser intensity of I = 5 x 10 14 
W/cm 2 . Sub-plots (a) and (d) present the dipole moment, (b) and (e) the dipole acceleration, and (c) and (f) the spectral density. There is 
remarkably close agreement between the TDDFT and SAE models. The dipole (center-of-charge) motion is consistently reproduced in phase 
and amplitude as can be seen from the dipole moment sub-plots. The irregular fluctuations in the dipole acceleration are also remarkably 
similar with a very efficient conversion of the 7th harmonic, extending to the 13th harmonic at the higher intensity. The oscillations in (e) for 
the SAE model after the pulse are due to Is — 2p fluorescence. There is negligible ionization for both models at the two intensities considered: 
less than 0.1% for I = 1 x 10 14 W/cm 2 and less than 2% for / = 5 x 10 14 W/cm 2 . Instead of a classical plateau and cutoff, we have a dip and 
secondary maximum below the classical cut-off lu c : equation J28l >, 



The theory and applications of S-splines are well known l32l,l33tl . but we briefly summarize their use for this problem. B-spline 
functions are localized overlapping piecewise polynomials designed to reproduce the radial oscillations of the wavepacket. These 
overlaps give rise to a narrow-banded symmetric structure in S, and by design, these elements can be evaluated exactly by Gauss- 
Legendre quadrature. A suitable choice of spline order, k, is governed by the Hamiltonian, and the radial space, < r < r max , 
subdivided in sectors or scaled to vary the density of points, accordingly. Specifically, we discovered that k = 9 order functions 
combine the advantages of low bandwidth and accurate representation of the wavefunction oscillations. 

The ground-state eigenvector is calculated and normalized to provide the stationary initial state c(0). Its evolution is then 
calculated by numerical solution of the system of equations d20l i. This is not a trivial task given the scale of the problem (the 
size of the interaction matrix) and the instabilities of the equations. It is well known that the choice of gauge has a strong effect 
on the dynamical terms in the matrices. Not surprisingly, this computational challenge has attracted a great deal of attention. 
While there exists a variety of tried and tested strategies for such problems, we have developed our own computer codes that 
treat the problem in complementary ways for both length and velocity gauges. This double approach provides both a numerical 
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FIG. 2: (Color online) Time-frequency analysis of high-order harmonic emission in a helium atom exposed to a five cycle (6.4fs) linearly 
polarized laser pulse of wavelength A = 390nm for two different models: a TDDFT model (left-hand plots) and a SAE model (right-hand 
plots). The spectral density of the STFT fequationll26t) is plotted as a function of the auxiliary time r; high density is colored red and low 
density blue. Two laser intensities were considered: sub-plots (a) and (b) correspond to / = 1 X 10 14 W/cm 2 while sub-plots (c) and (d) 
correspond to / = 5 x 10 14 W/cm 2 . 



validation for the physical results, and allows us to select the more efficient method when required. 



All calculations were performed with the PYPROP package H34I1 . which is designed to provide a general interface for solution 
of the time-dependent Schrodinger equation. Discretization schemes and propagation methods are implemented following a 
standard interface, which enables users to test different methods without requiring detailed knowledge about the implementation. 
PYPROP is written in Python, while all the computationally-intensive routines have been developed in C++, using the blitz++ 
array library ll35ll . Depending on the problem at hand, the size of the grid and number of employed grid points will vary. For 
the range of angular momenta, we found that l max = 20 was more than sufficient for convergence in all cases considered here. 
In addition we found that 1050 B-splines were sufficient for convergence. Time propagation was carried out using the Arnoldi 
method JH. 
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FIG. 3: (Color online) The effect of carrier-envelope phase (CEP), <p, on harmonic generation in a helium atom exposed to a five cycle (6.4fs) 
linearly polarized laser pulse of wavelength A = 390nm and intensity / = 1 x 10 14 W/cm 2 for two different models: a TDDFT model (left-hand 
plots) and a SAE model (right-hand plots). A pulse envelope defined by equation l llOt was used and two CEP values were considered: if = 
(black lines) and tp = tv/2 (red lines). Sub-plots (a) and (d) present the dipole moment, (b) and (e) the dipole acceleration, and (c) and (f) the 
spectral density. We note that, for both models, there is very little effect of CEP in the time-integrated spectral density. 



IV. FREQUENCY ANALYSIS 



The dipole moment is defined as 



d(t) = - J n(r,t)zd 3 r, 

= -2{m\Am)- 



(23) 



According to Larmor's formula, the radiated power is proportional to the square of the dipole acceleration. This can be obtained 
using Ehrenfest's theorem as ll37ll 



d(t) = -2(V>(i)l [H, [H,z]]m))- 



The spectral density is then obtained as follows 



S(u) = 



d{t)e lult dt 



(24) 



(25) 



8 




Harmonic order Harmonic order 
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FIG. 4: (Color online) The effect of carrier-envelope phase (CEP), <p, on harmonic generation in a helium atom exposed to a five cycle (6.4fs) 
linearly polarized laser pulse of wavelength A = 390nm and intensity 1 = 5 X 10 14 W/cm 2 for two different models: a TDDFT model (left-hand 
plots) and a SAE model (right-hand plots). A pulse envelope defined by equation l llOt was used and two CEP values were considered: if = 
(black lines) and tp = tv/2 (red lines). Sub-plots (a) and (d) present the dipole moment, (b) and (e) the dipole acceleration, and (c) and (f) the 
spectral density. The spectral density obtained from the TDDFT model exhibits a strong dependence on carrier phase, both below and above 
the secondary maximum. This behavior is consistent with a fully-correlated numerical simulation for helium in the infrared |6], Conversely, 
the spectral density obtained from the SAE model lacks this feature. 



To analyze the frequency content of a signal over time, both time and frequency information is needed simultaneously, i.e. a 
time-frequency analysis provides valuable insight. Here we use the Short-Time Fourier Transform (STFT), which has been 
found to recover dominant frequencies within a signal with reliable accuracy 13811 . The spectral density of the STFT is given 
by S 



F(lu,t) = 



+oo 



d\t) h(t 



'dt 



(26) 



where u> and r represent the analyzing frequency and time-shift respectively, and we chose to employ a Hann window function 
given by lf40h 



h(t) 



1 



2irt 
Tw 



(27) 



where Tw is the window length. In our case, given the timing of the radiation bursts, a natural choice for Tw was the optical 
period. 
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V. RESULTS AND DISCUSSION 

We present results for the dipole moment, dipole acceleration and corresponding power spectra for the two models at two dif- 
ferent laser intensities in FigQ] In all simulations a laser wavelength of A = 390nm was considered, corresponding to frequency- 
doubled Ti:sapphire wavelength. The two laser intensities considered are / = 1 x 10 14 and 5 x 10 14 W/cm 2 . For the TDDFT 
method, the ground-state energy, calculated by propagating in imaginary time, was found to be 2.7240 a.u. The ground-state of 
singly ionized helium is 1.8068 a.u., and thus the ionisation energy in the TDDFT model is I p = 0.9172 a.u. In the SAE model 
I p = 0.9034 a.u. which coincides with the ionization energy of the TDDFT model. The good agreement is perhaps surprising 
given the very different forms of the Hamiltonian and their respective treatment of correlation. We also calculated the static 
polarizability a, which describes the quadratic (lowest-order non-zero) Stark shift. For the TDDFT model presented we find 
a = 1.76, which is in good agreement with experiment a = 1.38 and other TDDFT calculations lf27tl . 

Let us now consider the center-of-charge oscillations in Figs, da) andOJd). At both intensities, the motion is dominated by 
the fundamental dipole motion. This manifests itself in the agreement of the fundamental frequency in the harmonic spectra for 
both models. This collective motion, in a many-body system, is expected of a tightly-bound oscillator under a low-frequency 
external field (A = 390nm corresponds to u>l — 0.1168 a.u.). In general, the center-of-charge motion is weakly dependent on 
the internal forces, and would not be expected to show dynamic correlation effects. Thus the different models of correlation in 
TDDFT and SAE models are not significant when the fundamental mode is observed. 

Deviations between TDDFT and SAE are apparent with high-frequency quivering at the peak of the external field in Figs.[TJa) 
and [TJd), and are revealed when the dipole mode is removed. The small amplitude oscillations are most pronounced at the 
peak of the electric field and occur near the nucleus where the restoring force is strongest and the electron acceleration is most 
intense. The same trend carries over to the dipole accelerations in Figs.[Tfb) andQJd), resulting in a more intense spectral density 
as shown in Figs. [He) andQJf)- 

Three features typical of HHG are J4ll I42L l43fl : (i) low order harmonics, which arise from transitions from bound excited 
states to the ground-state, show a rapid decrease in intensity as expected from a perturbative process; (ii) a (short) plateau region, 
caused by transitions from the continuum, of relatively constant intensities; and (iii) an abrupt cut-off, at a frequency 

lu c ~ Ip + 3.2U P , (28) 

where Ip — 0.9040 a.u. is the ionization potential. Rather than a plateau, we observe a dip followed by a secondary maximum. 
As the laser intensity is increased, this secondary maximum extends to higher frequency and broadens. Comparison with the 
classical cut-off formula, equation (|28]), for / = 1 x 10 14 and I = 5x 10 14 W/cm 2 , gives oj c ~ 9ujl and lj c ~ 15u>l, respectively. 
Fig.Q]illustrates, that both the TDDFT and SAE calculations falls short of these harmonics, with secondary maxima nearer the 
7th and 1 1th harmonics. On further investigation, we found this attenuation is not a modulation feature due to the pulse shape, 
in which the intensity variation is rapidly varying. A much longer pulse (T — 25.6fs) gives similar results for the drop-off 
frequency. 

Given that TDDFT and SAE provide conflicting descriptions of the asymptotic effective potential, it is clear that long-range 
correlations in the two models will differ and the ionization rates, for example, will disagree. At both these laser intensities we 
found that ionization at the end of the pulse was negligible in the two models (< 2%). Moreover, in the case of short-range 
behavior, the correlation terms in each model will be negligible in comparison with the Coulomb singularity that gives rise to the 
Kato cusp in the wavefunction. Essentially correlation effects are dominated by the nuclear singularity in this region. According 
to equation d24l > the highly-singular acceleration operator is strongly peaked in this region and thus the effective Hamiltonians 
are single-electron (uncorrelated hydrogenic) expressions at the highest frequency. Thus the TDDFT and SAE models share the 
same essential features in the lowest and highest frequency emission. At low frequency we have a collective low-frequency mode 
of the orbital, while the highest frequency emission is essentially single-electron hydrogenic behavior. So, at the extremes of the 
spectrum this again dominates correlation effects, while in the intermediate region (Fig.[TJc), harmonics 3, 5 and 7) correlation 
corrections are significant. 

At much higher intensity, I = 5 x 10 14 W/cm 2 , the laser field, Eq w 0.12, competes with the binding potential and the 
correlation energy. The degree of atomic excitation is greatly increased and we would expect significant disparity between the 
models. The SAE potential is weaker than the equivalent TDDFT effective potential. Nonetheless, ionization is less than a few 
per cent and Fig. [TJd) indicates that the center-of-charge motion is more or less the same. The dipole acceleration in Fig.[TJe) 
is now smeared across the duration of the pulse, though still indicating prominent bursts at the half-cycle turning points of 
the moment. The time-integrated spectral density in Fig. |TJf) shows the matching of the fundamental in both models, as was 
observed at the lower intensity, and again a stronger yield of lower harmonics in the TDDFT model. The secondary maximum 
is as strong as the fundamental peak. For example with a five-fold increase in intensity, the 9th harmonic is over 100 times more 
intense, as shown in Figs.QJc) andQJf)- The long slow decay of the SAE spectrum in the drop-off is consistent with a lower 
ionization threshold for this model. 
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The time-integrated spectral densities do not reveal the attosecond dynamics of the correlation. In Fig. [2] we present an 
analysis of each burst of harmonics through the STFT. The 5 cycle pulse provides an interaction time of T = 6.4fs. At the lower 
intensity, Fig.|2a), the burst created at 2.25fs is correlated to the turning point in Fig. [Ha). The subsequent burst, half a cycle 
later at 2.9fs shows a bifurcation at the higher frequencies which coincides with the double-peak in the dipole acceleration in 
Fig. Q2b). The irregularities in Fig. Q~fe) at the higher intensity are reflected in fringes within the main bursts in Figs. |2jc) and 
I2d) above. The peak emission time is also broadened and slightly delayed. In Fig. EJd), the SAE results show an extended 
fluorescence signal due to the 2p — Is emission after the passage of the pulse. At the higher intensity, Fig.|2jc) we notice the 
enhancement in high-frequencies at the expense of the lower frequencies, in accordance with Fig.Qlc). However, the timing of 
the peak bursts in Fig. [2c) are slightly delayed compared with Fig. |2ja), for example, with the TDDFT model. This delay is not 
observed for the SAE model. 

Baltuska et al [!6i] found that HHG was very sensitive to the carrier-envelope phase. Their experiments with helium, using 
the fundamental TkSapphire mode in the infrared A = 800nm, showed that, provided the CEP was chosen to coincide with the 
peak of the envelope, the spectrum generated in the cutoff region loses the odd-order harmonics. These results were supported 
by theoretical simulations in the same letter. We examined whether this holds true for our two models at the frequency-doubled 
laser wavelength: A = 390nm. The results are presented in Figs.|3]and|4] Firstly, referring to Fig. [3] the time-integrated spectral 
density, subfigures (c) and (f), show very little sensitivity to the phase ip. 

At the higher intensity the sensitivity to CEP is apparent in the TDDFT method but not the SAE results. The TDDFT 
simulations, Fig. Hfc), are consistent with experiment and simulations (fully-dimensional, fully-correlated) reported for A ~ 
800nm [6]. Beyond the cutoff it is clearly evident that the spectrum is continuous when (p = ir/2 and harmonic when ip = 0. 
On the other hand, the SAE simulations do not exhibit this behavior. Furthermore, the highest harmonic produced by the SAE 
model at I = 1 x 10 14 W/cm 2 is 13th-order and at 7 = 5 x 10 14 W/cm 2 is 23rd-order, while the TDDFT results extend well 
beyond this range. On this basis, it suggests that, for these parameters, a more complete model of correlation is provided by 
the TDDFT model. Of course, confirmation of this finding from a full-dimensional two-electron simulation and experimental 
measurements would be important. 



VI. CONCLUSIONS 



In this paper, we have investigated HHG in helium subjected to short intense laser pulses, using two models of correlation: 
a TDDFT approach and a SAE model. We presented harmonic yields for two laser intensities (7 = 1 x 10 14 W/cm 2 and I = 
5 x 10 14 W/cm 2 ), for the frequency-doubled TkSapphire laser: A = 390nm. We found that the linear response in both the 
TDDFT and SAE models are in good agreement, despite differences in the asymptotic behavior of their effective potential. This 
is reflected in the fundamental frequency in the harmonic spectra. We observed a secondary maxima below the classical cutoff 
frequency and this was found to be independent of the pulse duration. We investigated the effect of changing the CEP of the laser 
pulse on the harmonic yield, and found that the sensitivity of the harmonic spectra to the CEP is apparent at the higher intensity 
for the TDDFT calculation, which is consistent with previous experimental results. Finally, we presented a time-frequency 
representation illustrating the instants of high-frequency ultra-short bursts of light. 

Further work will involve optimization of both the TDDFT and SAE effective potentials to give a more accurate and realistic 
description of the system under investigation. It would be of interest to study pulse envelope effects on the HHG process using 
these models. 
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